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Abstract. A thorough analysis of the OGLE-II time-series photometry of the Large Magellanic Cloud bar sup- 
plemented by similar data from the MACHO database led us to the discovery of three /3 Cephei-type stars. 
These are the first known extragalactic f3 Cephei-type stars. Two of the three stars are multiperiodic. Two stars 
have inferred masses of about 10 Mq while the third is about 2 mag brighter and at least twice as massive. All 
three variables are located in or very close to the massive and young LMC associations (LH41, 59 and 81). It is 
therefore very probable that the variables have higher than average metallicities. This would reconcile our finding 
with theoretical predictions of the shape and location of the /3 Cephei instability strip in the H-R diagram. The 
low number of /3 Cephei stars found in the LMC is another observational confirmation of strong dependence of 
the mechanism driving pulsations in these variables on metallicity. Follow-up spectroscopic determination of the 
metallicities in the discovered variables will provide a good test for the theory of pulsational stability in massive 
main-sequence stars. 
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1. Introduction 

Since the early 90-ties it is known that pulsations in 
P Cephei stars are driven by k mechanism operating 
at a temperature around 2x10^ K where the opacity 
bump occurs (Cox et al. 1992; Moskalik & Dziembowski 



1992; Kiriakidis et al. 1992; Dziembowski & Pamyatnykh 
1995; Gautschy & Saio 1993). The bump is caused by 



a large number of bound-bound transitions in the iron- 
group ions. As a consequence, the area of the instabil- 
ity strip in the Hertzsprung-Russell (H-R) diagram should 
strongly depend on metallicity (Moskalik & Dziembowski 
1992; Pamyatnykh 199£). The instability in models of 



f3 Cephei stars results from a small domination of driving 
over damping (see, e.g., Moskalik & Dziembowski 1992D , 
thus the subtle differences in both the input physics and 
the modeling may lead to considerable differences in the 
predicted position of the instability strip. For instance, 
for solar metallicity, Pamyatnykh ( p_999 ) finds that the 
f3 Cephei instability strip (hereafter BCIS) has no upper 
edge for the most massive stars and extends beyond the 
core-hydrogen burning phase of evolution. Moreover, for 
the heavy elements abundance parameter Z = 0.01, his 
BCIS shrinks down to a small region at high luminosities, 
where /3 Cephei stars have not been found. On the other 



the much narrower BCIS with an upper boundary for M « 
20 Mq . In their calculations the BCIS disappears at Z as 
low as 0.005. 

In this context, the observational study of f3 Cephei 
star pulsations in objects of different metallicity is of great 
importance. I t was already pointed out by Sterken & 
Jerzykiewicz (1988) that with their lower than Galactic 
metallicities and relatively small interstellar absorption, 
the Large and Small Magellanic Clouds (LMC and SMC) 
are among the best objects for such a study. Although typ- 
ical photometric amplitudes of f3 Cephei stars are of the 
order of 0.01 mag and these stars should not be brighter 
than ^^14 mag in the LMC and even fainter in the SMC, 
with the modern techniques allowing the CCD photometry 
in crowded fields, the detection of these stars is certainly 
within reach. 

Searches for (3 Cephei stars in the Magellanic Clouds 
had already been undertaken. The first was carried out by 



hand, for solar abundance, Deng & Xiong (2001) predict 
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Sterken & Jerzykiewicz (1988). By means of the photoelec- 
tric photometry these authors studied six late O / early B- 
type stars. For the same purpose, Kubiak (1990) searched 
the young LMC cluster NGC 1712. These searches resulted 
in the discovery of some variables, but none of a convinc- 
ing (3 Cephei-type pulsation mainly because of the small 
statistical sample observed. Ano ther search for f] Cephei 
stars was perfor med by Balona (1992; 1993 ) and Balona 
& Jerzykiewicz ( |l993D in NGC 2004 and 2100 in the LMC 
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and in NGC330 in the SMC. A simil ar CC D search was 
also carried out by Kjeldsen & Baadc ( |1994D in NGC 2122 
in the LMC and NGC 371 in the SMC. No variable of 
(3 Cephei- type was found by these authors. 

The reduction of the OGLE-II data (Udalski et 
al. 1997) for the Magellanic Clouds by means of the 



Difference Image Analysis (hereafter DIA, Wozniak ^000| ; 



et al. 



2001a) yielded the photometry for over 
the LMC 



Zebruh 

53,000 candidate variables in the LMC and 15,000 in 
the SMC. A catalog including these stars was recently 
made available to the astronomical community (Zebruh et 
al. |200lH) . EarUer, the OGLE BVIq photometry of about 
6.7 X 10^ stars in the LMC (Udalski et al. |2000D and 2.2 
X 10^ stars in the SMC (Udalski et al. |199^ ) was pub- 



lished. At the same time, the even larger database of the 
MACHO observations of 7.3 x lO"^ stars in both Clouds 
and in two filter s {V , R) has become available (AUsnian 
& Axelrod |200ll) . 

These databases offer an unprecedented opportunity to 
search for the presence of /3 Cephei-type stars in the LMC 
and SMC. Combined with the detailed metallicity deter- 
minations, this result will put strong constraints on the 
theoretical stability predictions. As we shall show later, 
the accuracy of the photometry is good enough to detect 
all /3 Cephei-type pulsations in the LMC stars with semi- 
amplitudes exceeding 3 to 5 mmag, depending on the star 
brightness. 

In this paper, intended to be the first of a series, we 
present the results of the search for (3 Cephei stars among 
the LMC bar stars using the OGLE-II DIA photometry. A 
similar analysis for the SMC stars will be published sepa- 
rately. The details of the analysis and the selection of stars 
is described in Sect. 2. Section 3 presents the main results 
of the search. The possible connection of the variables we 
found with the LMC clusters and associations is verified in 
Sect. 4. We also briefly discuss the variables detected dur- 
ing previous attempts to flnd fi Cephei stars in the LMC 
(Sect. 5). The consequences of our finding are discussed 
in detail in Sect. 6, while a short summary and our plans 
concerning the next papers of the series are given in the 
last section. 



2. Data selection and analysis 

The output of the DIA is the photometry of the vari- 
able stars only (strictly speaking the candidates for vari- 
able stars, as the DIA produces also artifacts), therefore 
it was an obvious choice to take the OGLE-II catalog of 
about 53,400 candidate variables in the LMC (Zebruh et 
al. |2001bD as the main source of the time-series photomet- 
ric data. For the details of the OGLE-II observations and 



the DIA reductions we direct the reader to the paper of 
Zebruh et al. ( 2001a ). We only mention that the OGLE- 
II observations of the LMC cover twenty-one 14f2 x 56f9 
fields in the bar of this galaxy. For a single star, about 
30, 40, and 400 observations are typically available in the 
B, V and Ic bands, respectively. The OGLE-II observa- 
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Fig. 1. Colour-magnitude diagram for '-^16,000 stars in 
the OGLE-II LMC field SCI, rando mly selected from the 
photometry of Udalski et al. (2000). The three (3 Cephei 



stars, VI, V2, and V3, are shown as filled circles, the other 
two short-period variables, V4 and V5, as open circles. For 
reference, the location of the luminosity class V B-type 
stars is shown considering an LMC distance modulus of 
18.5 mag and Ay = 0.4 mag. The inclined solid line is the 
reddening line for Ay/E{V - Iq) = 2.49 (Stanek 11991). 



The dashed line shows the limits of the search for /? Cephei 
variables. 



tions span at most 3.5 yr, from the beginning of 1997 to 
mid-2000. 



The MACHO data archive (AUsman & Axelrod pOOlj ) 
served as the second source for the time-series data. The 
data span about 8 years between mid-1992 and the be- 
ginning of 2000, and are available in two passbands: blue 
(440-590 nm) and red (590-780 nm). Henceforth we shall 
refer to these passbands as Vm and i?M, respectively. The 
MACHO observations cover a much larger area around the 
LMC than the OGLE-II ones do. However, the database is 
at present accessible in an unsuitable way for our search. 
We therefore use here the MACHO archive only for stars 
we found particularly interesting from the analysis of 
the OGLE-II data and for variables detected in previous 
searches (see Sect. 5). A thorough search for /? Cephei stars 
in the MACHO archive will be done as soon as enhanced 
search capabilities become available. 

Prior to the analysis, the observations with larger pho- 
tometric errors have been rejected from the MACHO data. 
In addition, some outliers were removed. Finally, the he- 
liocentric corrections were applied to the MACHO epochs. 
Since the published epoc hs correspond to the beginning of 
an exposure (Cook ^002 ), 150 seconds correction (half the 
exposure time) was added to them. The OGLE-II epochs 



A. Pigulski, Z. Kolaczkowski: /3 Cephei stars in the LMC bar 



3 



have been corrected for the 'drift-scan effect' according to 
the equation given by Zebrun et al. (2001b). 

The coldest and faintest Galactic f3 Cephei stars 
have MK type of about B2.5V, corresponding to My « 
—2. Taking into account the LMC distance modulus of 
18.5 mag, we get V — 16.5 mag. The average E{B — V) 
colour excess for the LMC is about 0.13 mag (Massey et 
al. 1995), yielding Ay « 0.4 mag. However, the visual ex- 
tinction is sometimes much larger, especially for hot stars 
(Zaritsky 1999). This is obvious, because there is usually a 
large amount of interstellar matter around young hot ob- 
jects. In order to take that into account, we decided to add 
an 1 mag margin to the limiting magnitude of our search. 
Considering also the effect of reddening on the photomet- 
ric indices, we have finally chosen stars with y < 18 mag 
and {V — Ic) < 0.5 mag (see Fig. |l|) to search for the pres- 
ence of (3 Cephei-type variability. A total of 5204 stars in 
the LMC were extracted from the catalogue of Zebrun et 
al. ( ^OOlbl ). 

For all these stars we performed frequency analysis by 
means of the AoV periodogram of Schwarzenberg-Czerny 
(1996) in the range between and 20 d~^. The frequencies 
of the highest peaks in all periodograms were then sorted 
and the photometry of stars with the highest frequencies 
was examined visually. 

In the classical approach, (3 Cephei stars are recog- 
nized as early B-type stars with p eriods shorter than 
0.3 d (Sterken & Jerzykiewicz 1993 ). Since it seems that 
(3 Cephei stars with period s slig htly longer than 0.3 d ex- 
ist (Krzesinski & Pigulski 1997| ), we decided to consider 
all stars showing periods shorter than 0.35 d as potential 
variables of this type. 

As a result of our periodogram analysis, we found three 
stars which we regard to be the LMC (3 Cephei stars. Two 
other short-period objects were found among early-type 
variables. They could also be (3 Cephei stars, but because 
another interpretation of their variability is more likely, 
we describe them only briefly in Sect. 3.4. 

For the sake of simplicity we designate the (3 Cephei 
stars as VI to V3 and list them in Table 1. The cross- 
identifications with the OGLE-II and MACHO designa- 
tions are also given. The finding charts shown in Fig. ^ 
were taken from the OGLE-II database. 



3. (3 Cephei stars 

3.1. VI = OGLE 053446.82- 694209.8 
MACHO 81.8881.161 



■ - * *A ^ 



N 




Fig. 2. OGLE-II maps (lfl7 x lfl7) for the three 
l3 Cephei stars (V1-V3) in the LMC bar. The J2000 coor- 
dinates of the stars are coded within the OGLE-II names 
(Table 1). 
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Fig. 3. Fourier periodograms of the OGLE-II 1997-2000 
/c-filter observations of VI. The upper panel shows the 
periodogram of the original data; the lower one, after 
prewhitening with / — 4.05160 d^^. 



The star is located in the field SCI and 81 of the OGLE-II 
and MACHO, respectively. Fourier analysis of the OGLE- 
II Ic observations of this star (Fig. ^) revealed a single 
frequency / = 4.05160 d~^. No other periodicity above 
the noise level (~5 mmag) has been found. Practically 
the same result, that is, a monoperiodic signal, has been 
detected in the MACHO data except that after prewhiten- 
ing with /, a small peak appeared at frequency of 1 (side- 
real day)~^. It is presumably of instrumental origin. The 



star has virtually the same amplitudes in the Vm and Rm 
bands (see Table 2). Since the MACHO data span a wider 
time interval than the OGLE-II ones do, the refined fre- 
quency, / = 4.051593 d~^, has been obtained by means of 
a non-linear least squares fitting of the MACHO data. The 
semi-amplitudes and phases of both datasets are given in 
Table 2. 
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Table 1. The (3 Cephei stars in the LMC bar 





OGLE-II 


OGLE-II 


MACHO 


V 


B-V 


V -Ic 


Name 


field 


name 


name 


[mag] 


[mag] 


[mag] 


VI 


LMC-SCl 


OGLE053446. 82-694209. 8 


81.8881.161 


16.691 


-0.109 


-0.089 


V2 


LMC-SC3 


OGLE052809. 21-694432.1 


77.7792.493 


16.748 


+0.003 


-0.103 


V3 


LMC-SC7 


OGLE051841.98-691051.9 




14.327 


-0.175 


-0.181 



3.2. V2 = OGLE 052809. 21- 
MACHO 77.7792.493 



694432.1 



The second /3 Cephei-type variable has been found in the 
OGLE-II field SC3 and in the MACHO field 77. It has 
almost the same colour and magnitude as VI (see Fig. ^ 
and Table 1). The Fourier periodogram of the OGLE-II 
Jc-filter data (Fig. ^) shows the strongest signal at fre- 
quency /i = 3.49767 d~^. After removing this signal the 
highest peak appears at a frequency fsd = 1.00333 d~^, 
which is close to 1 (sidereal day)~^. Usually, such a sig- 
nal appears in the periodogram when long-term irregu- 
lar changes are superimposed. Since it appears for V2 in 
all three datasets, it can be caused by a low-frequency 
intrinsic variation. In fact, changes in the mean magni- 
tude of this star are well seen in the MACHO data: dur- 
ing the first three years of observations the mean bright- 
ness of V2 faded by about 0.02 mag. After removing /i 
and the low-frequency signals, a second periodicity, /2 = 
3.68377 d~^, with a very small amplitude appears slightly 
above the noise level (Fig. ^). In order to make sure this 
is a real frequency, we calculated similar periodograms for 
the MACHO data (Fig. |b). Although barely visible in 
the Vm observations, /2 appears clearly in the Rm data. 
We conclude that the frequency is real, and thus V2 is a 
biperiodic variable. However, because of the low amplitude 
of this periodic signal, a few close peaks of similar height 
occur around /2 . The highest peak does not have the same 
frequency in all periodograms; it is equal to 3.6838, 3.6853, 
and 3.6847 d"! for the OGLE-II Iq, MACHO Vm and 
MACHO Rm data, respectively. In the final fit we used 
the frequency corresponding to the highest detected peak 
among the different periodograms (see Table 2), but we 
cannot be certain that our choice is correct. 



3.3. V3 = OGLE 051841.98- 691051.9 

The third star found to be an LMC (3 Cephei star is about 
2 mag brighter than VI and V2 (Fig. |l|) and shows a very 
complicated frequency pattern. As many as five frequen- 
cies with signal-to-noise (S/N) ratio larger than 4 have 
been detected (see Sect. 3.5 for the definition of S/N). 
Figure ^ shows consecutive steps of prewhitening for the 
OGLE-II Iq data. One can see that all terms have very 
low amplitudes (Table 2). The last detected frequency is 
rather doubtful. Moreover, the frequency 2.0055 d^^ is 
very close to 2 (sidereal days)^^ which often appears in 
periodograms of bright stars. It is thus probably an ar- 
tifact. However, the other low- frequency term, with P « 




2 3 4 5 6 

FREQUENCY [d"'] 



3 9 10 



Fig. 4. a. Fourier periodograms of the OGLE-II 1997- 
2000 /c-filter observations of V2. From top to bottom 
the panels show the periodograms (i) of the original data, 
(ii) after prewhitening with /i = 3.49767 d~^, (iii) after 
prewhitening with /i and fsd — 1.00333 d~^, (iv) after 
removing /i, fsd, and /2 = 3.68377 d^^. b. Fourier pe- 
riodograms of the MACHO Vm (top) and Rm (bottom) 
data after prewhitening with the low-frequency variations 
and /i. The triangles indicate the location of /2 and its 
1 d^^ aliases. 



0.5938 d, can be real. According to the recent calculations 
of Pamyatnykh ( |l999| ) , massive stars in this region of the 
H-R diagram have unstable modes of such long periods, 
thus the low-frequency variation may also represent pul- 
sation. 
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Table 2. The results of the sme curve fits to the photometric data of V1-V3. iVobs is the number of individual 
observations and RSD is the residual standard deviation. See Sect. 3.5 for the definition of the signal-to-noise (S/N) 
ratio. 





Frequency 


Period 


Data 






Semi-amplitude 


HJD of the time 


RSD 




Star 




[d] 


source 


Filter 


iVobs 


[mmag] 


of maximum light^ 


[mmag] 


S/N 


VI 


4.051593 


0.2468165 


OGLE-II 


B 


22 


45.1 ± 4.8 


51 137.1511 ± .0043 


15.5 


7.8 








OGLE-II 


V 


45 


37.7 ± 3.1 


51 058.4247 ± .0036 


15.5 


9.3 








OGLE-II 


Ic 


350 


38.9 ± 1.0 


51027.8218 ± .0010 


12.7 


32.3 








MACHO 


Vm 


816 


39.7 ± 0.9 


49 947.2654 ± .0009 


18.1 


23.3 








MACHO 


Rm 


767 


37.2 ± 0.9 


49 896.1763 ± .0010 


17.9 


32.7 



V2 3.497909 0.2858851 OGLE-II Ic 501 12.4 ± 1.1 50 837.3568 ± .0041 17.5 10.4 

MACHO Vm 1351 15.2 ± 0.8 49 938.5438 ± .0021 18.9 16.3 

MACHO Rm 1188 13.6 ± 0.8 49 915.6678 ± .0027 19.9 13.0 



3.685343 0.2713452 OGLE-II Ic 501 4.0 ± 1.1 50 837.3874 ± .0118 17.5 3.3 

MACHO Vm 1351 3.6 ± 0.7 49 938.6814 ± .0089 18.9 3.9 

MACHO Rm 1188 3.9 ± 0.8 49 915.6335 ± .0093 19.9 3.7 



V3 5.179046 0.1930858 OGLE-II Ic 471 5.1 ± 0.5 50 852.1238 ± .0027 6.8 9.3 



3.495009 


0.2861223 


OGLE-II 


Ic 


471 


3.9 ± 0.4 


50 852.0745 ± 


.0054 


6.8 


7.1 


2.005502 


0.4986283 


OGLE-II 


Ic 


471 


4.5 ± 0.5 


50 852.2471 ± 


.0092 


6.8 


8.2 


1.684015 


0.5938189 


OGLE-II 


Ic 


471 


3.1 ± 0.4 


50 852.1539 ± 


.0137 


6.8 


5.6 


3.816132 


0.2620454 


OGLE-II 


Ic 


471 


2.3 ± 0.5 


50 851.9626 ± 


.0081 


6.8 


4.1 



^ The time is calculated to be as close as possible to the mean epoch of all observations of a given dataset. The first two digits 
of HJD, that is '24', have been subtracted. 



Unfortunately, the star was not observed by the 
MACHO and we cannot verify the reality of the modes 
independently, as we did for VI and V2. Also, there are 
not enough datapoints in the OGLE-II B and V band 
observations to detect reliably a similar variation. 



3.4. Other short-period variables 

As mentioned in Sect. 2, two other short-period variables 
may also be (3 Cephei-type stars. The first one, V4 — 
OGLE 051652.03-691405.2 = MACHO 79.5863.500 is lo- 
cated in the OGLE-II SC8 and MACHO 79 field. Its light 
curve is almost strictly sinusoidal in shape, has a period P 
= 0.260083 d, and the semi-amplitude of about 0.10 mag 
in all three OGLE-II filters. There are two characteristics 
which are difficult to explain if we assume the star to be a 
P Cephei variable. First, with V = 17.683, B-V = +0.066, 
and V—Ic = —0.008, the star falls in the colour-magnitude 
diagram about 1 mag below VI and V2 (Fig. g). Provided 
that V4 is not ~20 kpc behind the LMC — which is rather 
unlikely — this position corresponds to a mid-B-type star 
(Fig. |l|), too cool to have low-degree p modes excited. 
Second, for such a large amplitude, we would expect to 
see non-sinusoidal variations. For the two (3 Cephei stars 
with the largest amplitudes known, BW Vul and a Sco, 
the light-curve shows the so-called 'stillstand', a bu mp on 
the rising branch of th e light curve (Sterken et al 



Jerzykiewicz & Sterken 1984 ). Since this is not observed, 
classifying V4 as a Cephei variable is open to doubt. 



1986 



There is, however, an alternative explanation for the 
variability of this star. It could be a very rare contact 
or nearly contact binary consisting of two similar mid-B- 
type stars and seen at low inclination. Simple calculations 
indicate that with 2P « 0.52 d we can have a contact 
binary containing two main-sequence stars with masses of 
5-6 Mq. Such an object would have the position of V4 in 
Fig. |1|. 

The next star we suspected to be a /3 Cephei 
variable, V5 = OGLE 053811.81-704245.7 = MACHO 
11.9350.116, is very peculiar. It is found in the OGLE-II 
SC17 and MACHO 11 field. The star is fainter than VI 
and V2 {V — 17.509), and has much redder colours {B — V 
— -1-0.232, V — Ic = -1-0.351). However, it is located very 
close to the LMC nebula and cluster NGC 2075 = N213 
and in addition the SC17 field has a large reddening spread 
(see the colour-magnitude diagram presented by Udalski 
et al. 200C). Consequently, its colours could be accounted 
for by the reddening effect which will move the star, along 
the reddening line, close to VI and V2. 

If the star were a (3 Cephei variable, it would be indeed 
an unusual object. While the two detected frequencies (/i 
= 3.780 d-\ /2 = 3.777 d'^) are typical for a (3 Cephei 
star, the amplitudes are enormous. The amplitude of the 
/i term in B reaches 0.43 mag, exceeding about twice the 
amplitude in BW Vul. There is also a strong dependence 
of the amplitude on wavelength: in the Iq the amplitude is 
about half that in B. The second mode has smaller ampli- 
tudes, but still rather large for a f3 Cephei star, 0.24 mag 



6 



A. Pigulski, Z. Kolaczkowski: 13 Cephei stars in the LMC bar 




2^4567 

FREQUENCY [d"'] 

Fig. 5. Fourier periodograms of the OGLE-II 1997-2000 
/c-filter observations of V3. From top to bottom the pan- 
els show consecutive steps of prewhitening with frequen- 
cies (given in Table 2) indicated by triangles. 



1993 1994 1995 1996 1997 19 



2000 YEAR 



MACHO v.. 



ySS>" OGLE-II \r 



49000 49500 50000 50500 51000 

HJD 2400000. + 



Fig. 6. The OGLE-II Ic and MACHO Vm and Rm light 
curves of V5. Note the beating due to the presence of two 
periodicities with very close frequencies. Ordinate ticks 
are separated by 0.2 mag. 



in B and 0.10 mag in 7c- The two frequencies are very 
close to each other and the beating with a period longer 
than 300 days is clearly seen in the data (Fig. ^). We also 
found that both periods increase. The rate of change of P2 
= amounts to about 1.7 x 10~® and is seven times 
larger than that of Pi . As a result of the changes of the 
periods, the beat period shortened from 380 d in 1993 to 
about 317 d in 2000. 

Recalling the arguments given above for V4, we see 
that almost sinusoidal variations observed in V5 with even 
larger amplitudes than in V4 are difficult to reconcile with 



(3 Cephei-type variability. There is, however, another pos- 
sibility. The star can be an RR Lyrae variable of Bailey 
type c. The colours, periods, amplitudes and the shape 
of the light curve agree well with this type of variability. 
Moreover, RRc stars with very close periods have been re- 
cently found (Olech et al. [19994 |1999b| ; Moskah k|200q ). In 
the new classification scheme of Alcock et al. (200C), the 
star would be of RRl-wl type. However, V5 would then 
be too bright to be an LMC star. In accordance with its V 
magnitude, the star would be a foreground object located 
approximately half way to the LMC. There is another, in- 
direct argument in favour of this explanation: RR Lyraes, 
especially of the RRc type, are known to e xhibit large er- 
ratic period changes (see, e.g., Jurcsik et al. 2001 , Kopacki 
^OOlD which a re rat her not observed in the P Cephei stars 
(Jcrzykiewicz 1999 ). In any case, the star is very unusual 
and deserves further attention. 

As far as V5 is concerned, the low-resolution spec- 
troscopy would be sufficient to distinguish between the two 
possibilities. Even accurate Johnson U photometry, which 
is probably soon going to become available (Zaritsky et 
al. 1997) would be conclusive in this case. For V4 we need 
radial velocities to decide on the final classification. 



3.5. Detection threshold 

In view of the results presented here and the discussion 
that follows in Sect. 6, it is useful to estimate the detection 
threshold of the search. We therefore analyzed the light 
curves of 36 variables randomly selected from the cata- 



logue of Zebruh et al. (gOOlbj) with 14.5 <V < 19.1 mag. 
First, all periodic signals were removed from the data. 
Next, Fourier periodograms of the residuals were calcu- 
lated giving mean semi-amplitudes between and 10 d^^. 
We define this value as the noise (N). It is easy to find 
that if we define signal (S) as the height of the maximum 
peak in the Fourier periodogram, the signal-to-noise ratio 
(S/N) for a pure Gaussian noise would be (S/N)noisc ~ 3.0. 
Therefore, we have arbitrarily set our detection threshold 
to Dth = 4N. The thresholds defined in this way are shown 
in Fig. for the stars mentioned above. It can be seen from 
the figure that in the range of magnitudes where the LMC 
(3 Cephei stars are expected to occur, we should be able to 
detect amplitudes as low as 3-5 mmag. This agrees very 
well with the results given in Sect. 3.1-3.3 and Table 2. 
Since the OGLE-II DIA photometry practically does not 
suffer from crowding, the relation shown in Fig. ^ should 
apply to almost all variables. 

However, we have to remember that the analysis was 
carried out only on those stars, already found as variables 
by the DIA. In order to check whether the DIA missed 
some small-amplitude variables which could be detected 
by periodogram analysis, we calculated S/N for all can- 
didate variables in the SCI field using their Fourier pe- 
riodograms in the range between and 10 d^^. Prior to 
analysis, we additionally filtered the observations using a 
4(j clipping in order to avoid the infiuence of outliers. The 
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Several compact groups within LH81 are usually distin- 
guished as open clusters; the brightest one is the OGLE- 
CL-L MC60 5 (Pietrzyiiski et al. |l999| ) = BCDSP 8 (Bica 



16 17 
V [mag] 

Fig. 7. The detection threshold, Dth, for the OGLE-II 
DIA /c-filter photometry of 36 stars (dots) in the LMC 
as a function of the V magnitude. The continuous line is 
the polynomial fit to the points, the dashed line shows the 
detection threshold of 5 mmag. See text for the definition 
of Ah- 

S /N values calculated in this way have a very well defined 
limit of (S/N) w 3.6, independent of the magnitude. This 
shows that the DIA did not miss low-amplitude variables 
and its detection threshold is practically the same as that 
shown in Fig. |^. 

It is also interesting to compare the accuracy of the 
OGLE-II and MACHO data. As can be seen from the 
comparison of the RSD values in Table 2, the accuracy of 
a single observation is slightly better in the OGLE-II data. 
However, because there are usually more observations for 
a given star in the MACHO database, the noise levels in 
the periodograms are similar for both data sets. 

4. Relation with clusters and associations 

The P Cephei-type variability is confined to stars in the 
core-hydrogen burning phase of evolution or shortly be- 
yond this phase having masses larger than ^7Mq (see, 
e.g., Pamyatnykh 1999| ). From the comparison with evo- 
lutionary models it can easily be shown that a f3 Cephei 
star cannot be older than ~35 Myr. This is why almost all 
Galactic /3 Cephei stars can be identified as the members 
of the young open clusters or OB associations. We there- 
fore expect that the LMC /? Cephei-type stars can also be 
members of clusters or associations. 



4.1. VI and the LH 81 association 

Our first variable, VI, is located in the middle of the 
ring-shaped superbubble N 154 (Henize |1956|) = DEM 246 



(Davies et al. 1976) which is also a diffuse X-ray source 
(Dunne et al. |200lD . The two brightest parts of N 154, 
south-western BSDL2434 (Bica et al. |l999|) and the 
much brighter north-eastern LH 87n, with a bright knot 
NGC 2048, are well seen in Fig. ||. N 154 encomp asses two 
large OB associations: LH 81 (L ucke & Hodge 197C) ~ 
SL589 (Shapley & Lindsay |1963D and LH87 (see Fig. i). 



et al. 1996). The LH81 association has been studied by 
Massey et al. (2000) and is known to contain many very 
hot and massive objects including three Wolf-Rayet stars 
and two stars of spectral type 03-4. Figure |^ shows the 
colour-magnitude diagram for a region in LH81 with 
a radius of 2'.7 centered at (a,(5)2ooo.o (5^34'"39", 
— 69°43'11"). We see that the main sequence extends to 
stars as bright as eleventh magnitude and that VI is a 
likely member of the association. If this is true and the 
star formation was coeval, VI should not be olde r tha n 
5 Myr, the age of LH81 deduced by Massey et al. ( ^OOOl) . 



4.2. V2 and the LH 59 association 

V2 is located about 7' north-east from the center of 
the LH59 association (Fig. ||). The nearest clusters are 
the OGLE-CL-LMC499 = BSDL1861 and OGLE-CL- 
LMC505 = BSDL 1892. The colour-magnitude dia grams 
of the clusters provided by Pietrzynski et al. (1999) 



show that at least BSDL 1892 contains some hot stars. 
The ages of the two clusters estimated by Pietrzyiiski & 



Udalski ( pOOOD are 80 and 250 Myr for BSDL 1861 and 
1892, respectively. Thus, the clusters are much too old 
to contain (3 Cephei stars. However, LH59 seems to be 



slightly younger. Dieball & Grebel (200C) estimated the 



age of the three rich clusters embedded in this association 
(NGC 1969, 1971, 1972) to be in the range 40-70 Myr. 
Taking into account the errors in the age estimates and 
the possibility that the formation was not strictly coeval 
in such a large complex, the presence of a /3 Cephei star in 
this region might still be possible. Although we cannot be 
certain that V2 belongs to LH 59, the density of hot stars 
in this region is quite high. We conclude that V2 belongs 
either to the coronal region of LH 59 or to the LMC field, 
where, as fo und e.g. by Massey et al. (1995) and Holtzman 
et al. ( 1999| ), the star formation is still in progress. 

We also note that V2 has been detected as a far- 
ultraviolet source by the UIT telescop e. Its Bl-filter mag- 
nitude of 14.07 ± 0.07 (Parker et al. |l998| ) indicate that 
it is really an early B-type star. 

4.3. V3 and the LH 41 association 

It is, in turn, easy to indicate a connection of our third 
variable, V3, with an LMC association. The star is located 
in the rich LH41 association known to contain SDoradus, 
the extreme Luminous Blue Variable (LBV, Humphreys 
& Davidson 1994). The association (Fig. ^ contains also 
another probable LBV star, R85, a Wolf-Ra yet st ar, Brey 
21, a nd ma ny 0-type stars (Breysacher et al. 1999; Massey 
et al. 2000). Several young clusters including the compact 
SL360 = HD 35342 and the larger NGC 1910 = SL371, 
form the well-seen groups within the association. The 
whole region is surrounded by the N 119 — DEM 132 neb- 
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Fig. 8. Three 20' x 20' fragments of the red-filter POSS-II plates centered at VI (left), V2 (middle), and V3 (right). 
Associations, clusters and stars discussed in the text are labeled. White crosses mark the positions of /3 Cephei variables. 
Below, for each variable (encircled) the location of nearby bright, hot stars with < 18 mag and {V — Ic) < 0.5 mag, 
taken from the catalogue of Udalski et al. ( pOOO ), is shown on the same scale. North is up, east to the left. 
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Fig. 9. Colour-magnitude diagrams for the brightest stars of the LH 81 (left) and LH 41 (right) association. The 
photometry of stars shown as dots was taken from Udalski et al. (2000), those shown as crosses, from Massey et 
al. ( |2000| ). The Wolf-Rayet stars are shown as encircled open circles, LBVs, as encircled dots and f3 Cephei stars, as 
open circles, in addition to a label. Only stars with {B — V) < 1.0 mag and F < 18 mag are shown in the plots. 
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ula (Fig. H). There is no doubt that intensive star for- 
mation is taking place here. The age of the association 
estimated by Massey et al. (2000) is 4 Myr. The colour- 
magnitude diagram for the region of LH41 with a 4f5 
radius centered on (a,(5)2ooo.o = (5''18'"30", -69°13'15") 
is shown in Fig. |^. 

The UBV photometry of V3 has been obtained by 
Massey et al. ( |2000[ ). The star has been designated as 
LH41-64 by these authors and has V = 14.32, B - V ^ 
-0.16 and = -1.02 mag. Thus we get E{B-V) = 

0.14 mag and intrinsic colours: {B — V)o = —0.30 and 
{U — B)o = —1.12 mag which correspond rougly to a 
BO/09 type star. 

5. Notes on previous searches 

As mentioned in the Introduction, the previous searches 
for /3 Cephei stars in the LMC (Sterken & Jerzykiewicz 



1988 ; Kubiak 1990; Balona 1993 ) yielded some variables 



although none of them could be convincingly classified as 
a P Cephei variable. With the OGLE-II and MACHO pho- 
tometry at hand we can now try to verify their variability. 

Out of the six program and five comparison stars from 
Brunet et al. (1975) selected for observations by Sterken 
& Jerzykiewicz (1988), only star 90 and 134 are within 
the OGLE-II fields. Both stars were identified by the DIA 
as variables. However, the amplitudes are very small (less 
than 0.006 mag for star 90 and 0.01 mag for star 134) and 
are very likely of instrumental origin. On the other hand, 
the MACHO observed all these stars but 230. From our 
analysis of the ten stars we found that: (i) none of the stars 
shows short-period variability, (ii) stars 32 and 257 are 
variable, but the light changes are aperiodic and have the 
amplitudes below 0.2 mag, (iii) a 1-year periodicity is seen 
for four stars (90, 146, 230, and 134), but comparing the 
OGLE-II photometry of stars 90 and 134 we are certain 
that the variations arc of instrumental origin, (iv) stars 
269 and 271 are heavily blended. Concluding, none of the 
stars show P Cep hei-type variability. 

Kubiak ( 199C| ) found nine variables in the LMC cluster 
NGC 1712. Of these, seven (all but 5 and 8) are included 
in the MACHO field 44. Variable 6 is a blend and we 
could not cross-identify it properly, but the three nearby 
MACHO stars are constant. Of the remaining six stars, 
only variable 1 seems to be constant. Variable 2 shows a 
smooth increase and then decrease of brightness with an 
amplitude of 0.4 mag in Vm and 0.5 mag in Rm- Variables 
3 and 4 are eclipsing binaries with Algol- type light curves. 
The eclipses of equal depth (0.42 mag for variable 3 and 
0.17 mag for variable 4) indicate that the binaries have 
equal components. Variable 7 has a small amplitude and a 
probable period of 0.94 d, while variable 9 shows irregular 
variations mainly in the Vm band. 

Balona (1993) detected variability of seven stars in an- 
other LMC cluster, NGC 2004. Of those, five have pho- 
tometry in the MACHO database. We confirm the vari- 
ability of only two (Table 3). Star 87 is an eclipsing binary 
with a P Lyrae-type light curve and an orbital period of 



Table 3. Stars found to be variable in previous searches 
for P Cephei stars in the LMC identified with MACHO 
sources 

Star MACHO name Comment on variability 



Sterken & Jerzykiewicz ( 198J ) 



32 
90 
146 
210 
257 



18.2479.10 
79.4774.10 
49.6736.14 

8.8784.291 
50.9639.7 



Irregular variations 
Constant 
Constant 
Constant 

Irregular variations 



Kubiak (1990) 



1 44.1626.20 Constant 

2 44.1626.35 Long-period variable 

3 44.1626.46 Eclipsing, P = 2.291265 d 

4 44.1626.46 Eclipsing, P = 1.636929 d 
7 44.1747.28 P = 0.9403 d ? 

9 44.1868.9 Irregular 



Balona (1995) 



87 61.8192.120 Eclipsing, P = 1.662846 d 

241 61.8192.154 Irregular variations 

297 61.8192.55 Constant 

491 61.8191.12 Constant 

682 61.8191.81 Constant 



1.662846 d while star 241 is an irregular variable. The pe- 
riod of 1.835 d given by Balona (1993) for the latter star 
is not confirmed. Stars 297, 491, and 682 were found to be 
constant. 



The other LMC cluster searched by Balona (1993), 
NGC 2100, was not observed by either the MACHO or 
the OGLE-II. 



6. Discussion 

P Cephei-type stars have been finally found in the LMC. 
These are the first extragalactic stars of this type so far 
discovered. Photometric indices, magnitudes, amplitudes 
and periods of three variables fit the P Cephei charac- 
teristics. However, the spectroscopic confirmation and/or 
the spectral type of these variables would be, of course, 
desirable. 

Since there are over 27,000 early B-type stars in the 
LMC area covered by the OGLE-II observations, and the 
variability detection threshold was rather low, we may 
conclude that for the LMC the BCIS has been mapped. 
Only three P Cephei-type variables have been found. Two 
important facts arise from this finding. The first is the 
fraction of P Cephei stars with respect to the whole pop- 
ulation of early B-type stars, the other is their position in 
the colour-magnitude diagram. 
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6.1. Fraction of (3 Cephei stars with respect to early 
B-type stars in the LMC 

Practically all Galactic (3 Cephei stars have My between 
—4.5 and —1.5 mag. We take this range as a reference 
and in order to avoid selection effects we restrict the com- 
parison to stars in open clusters. Moreover, to make the 
comparison reliable, we need to account for different de- 
tection thresholds. A detailed analysis of the fraction of 
(3 Cephei stars in Galactic open clusters is going to be 
published soon (Pigulski et al., in preparation), here we 
mention only the preliminary result presented by Pigulski 
et al. (2002). These authors pointed out that the frac- 



tion of (3 Cephei stars in Galactic southern clusters is few 
times larger than in the northern ones, and, if related to all 
main-sequence stars falling into the My range given above 
it amounts to 35 ± 7% for the southern and 5.3 ± 1.7% for 
the northern clusters. Since southern clusters are closer to 
the Galactic centre, this was explained as being a conse- 
quence of the Galactic metallicity gradient. At the LMC 
distance, —4.5 < Afy < —1.5 mag corresponds roughly to 
14.3 < V < 17.3 mag. There are 27,663 main-sequence 
stars in the OGLE-II fields falling in this range of V mag- 
nitude. (As a main-sequence star we mean here a star with 
{y — Ic) < 0.5 mag; since contamination by foreground 
stars is very low in this part of the colour-magnitude dia- 
gram, it does not affect our estimate.) The three (3 Cephei 
stars we found, constitute 3/27663 « 0.011 ± 0.006% of 
early B-type stars (we assume a Poisson statistics to de- 
rive errors). Even if we take into account the fact that for 
searches in the Galactic open clusters roughly twice as low 
threshold is achieved as that seen in Fig. 0, so that some 
small-amplitude variables could have been missed in the 
LMC, the striking difference in the /? Cephei fraction be- 
tween the Galaxy and the LMC is evident. In other words, 
as the average LMC metallicity is [M/H] « -0.4 (or Z « 
0.007), this is another observational confirmation of the 
fact that the driving mechanism strongly depends on the 
metallicity. 

6.2. Position of (3 Cephei-type stars in the 

colour-magnitude diagram: a comparison with 
theoretical predictions for LMC metallicity 

Should we, however, expect to observe (3 Cephei stars at 
aU in the LMC if Z = 0.007? The two recent theoretical 



predictions (Pamyatnykh 1999; Deng & Xiong 2001) give 
different answers to this question. Pamyatnykh ( 1999| , see 
his Fig. 11) for Z = 0.01 does not find unstable modes for 
stars with masses smaller than 25 AIq. On the other hand, 
Deng & Xiong ( pOOl ) predict instability for stars with 
masses '^IOA/q down to Z = 0.005. The following reasons 
coming from our LMC results favour the predictions of 
Pamyatnykh (|1999D: 



shown for three values of Z: 0.01, 0.015 and 0.02. We 
see in this figure that for Z = 0.015 the lower part 
of the BCIS covers about half of the main-sequence 
width. Going towards higher masses the BCIS becomes 
narrower and then widens again. For Z — 0.01 there is 
no lower part at all and instability begins for masses 
'^25 Mq. It seems that for a value of Z between 0.01 
and 0.015 the BCIS splits in two regions, one for masses 
around 10 Mq and one over 25 Mq. This does corre- 
spond to the situation observed in the LMC. 
— In order to obtain the agreement mentioned above, 
we need, however, to have either a higher metallic- 
ity in the LMC (3 Cephei stars or the splitting of the 
BCIS describ ed ab ove, but for Z = 0.007. Although 
Pamyatnykh ( 1999|) , comparing the location of the hot 
border of the BCIS based on OP and OPAL opac- 
ities, estimates that an error in Z of about 0.003- 
0.005 can be involved, we think that the first possi- 
bility is more probable. We know from Sect. 4 that 
at least two, and perhaps all three f3 Cephei stars in 
the LMC are located in regions of violent star for- 
mation. The recent studies of the age-metallicity rela- 
tion for the LMC (Cohen |1982| ; Pagel & Tautvaisiene 



1998) and direct metallicity determinations in young 

Lambert 



objects (Olszewski et al. 
Jasniewicz & Thevenin 



1991 ; Luck 
1994; RoUeston et al. 



1992 



1996 



First, we find /3 Cephei stars both among low mass (8- 
10 Mq, VI and V2) and high ma ss (25 -30 Mq, V3) 
stars. In Fig. 11 of Pamyatnykh (|l999|) , the BCIS is 



Luck et al. 1998; Korn et al. 2000) indicate that al- 
though the mean LMC metalicity is about [M/H] « 
—0.4, a large spread in metallicities of young objects 
is observed. It is even possible that some objects have 
solar or nearly solar metallicity. In fact, it would be 
naive to believe that dozen megayears of the LMC his- 
tory did not leave any metallicity inhomogeinities. It is 
therefore very likely that our variables have a metallic- 
ity higher than the LMC average. Nevertheless, spec- 
troscopic determination of abundances in these stars 
would be highly desirable. 
— Why do we think the predictions of Deng & Xiong 
(2001) are incorrect? First, these authors do not find 
the instability among massive stars, so it would be dif- 
ficult to explain the variability of V3 with a lower than 
Galactic metallicity. Next, with instabilities predicted 
down to Z = 0.005, we should have found many more 
/3 Cephei stars than we did. Finally, Deng & Xiong 
( ^OOlD used the first 1992 version of the OPAL opac- 
ity tables with an analytical approximation which — as 
these authors write — gives as much as 10% differences 
with respect to the tabulated values. In view of the 
subtle domination of driving over damping, this could 
influence considerably the final result. 



7. Summary and future work 

The most important results obtained in this paper can be 
summarized as follows: 



1. Three (3 Cephei stars have been found in the LMC bar. 
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2. Two of them were found to have V magnitude corre- 
sponding roughly to the LMC B2-B2.5 stars, while the 
third one (V3) is as bright as a BO- type star. 

3. The three variables constitute a very small fraction 
(less than 0.02%) of all early B-type stars falling within 
the range of magnitudes where the Galactic stars of 
this type are distributed. This is much less than for 
the Galactic clusters and can be explained by a strong 
dependence of the driving mechanism on metallicity. 

4. The location and number of P Cephei variables in the 
LMC is in reasonable agree ment w ith the theoretical 
predictions of Pamyatnykh ( 1999 ) provided that the 
metallicity of these stars is slightly higher than the 
LMC average for young population, namely, between 
Z = 0.01 and 0.015 ([M/H] between -0.3 and -0.1). 

5. All three variables were found to be in or very close to 
regions of intensive star formation, implying that they 
indeed can have enhanced metallicities. 

6. If the low-frequency variation in V3 is confirmed, this 
can be an interesting example of a luminous (3 Cephei 
star in which modes with longer than typical periods 
are excited. The star could be therefore a good target 
for asteroseismology. 

The main goal of this series of papers is to give a thor- 
ough observational picture of the variability among the 
LMC and SMC stars lying in the upper main sequence, 
where the k mechanism related to the metallicity bump is 
at work. For this purpose we shall use, like in the present 
paper, the OGLE-II and MACHO photometric databases. 
In particular, we are going to focus on the study of slowly 
pulsating B stars (SPB), which, although fainter than 
P Cephei variables, have typical amplitudes large enough 
to be detected in the LMC and even in the SMC. Next, 
we are going to search bright B and A-type stars for vari- 
ability with periods of the order of 1 day and longer. This 
would probably shed some light on the presence of the g 
modes in luminous early- type stars. Finally, the variabil- 
ity of Be stars and hot eclipsing binaries will be studied 
as well. 
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